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Traffic overflowing a first-choice trunk group can be approximated 
accurately by a simple renewal process called an interrupted Poisson 
process-a Poisson process which is alternately turned on for an ex- 
ponentially distributed time and then turned off for another {independent) 
exponentially distributed time. The approximation is obtained by matching 
either the first two or three moments of an interrupted Poisson process to 
those of an overflow process. Numerical investigation of errors in the 
approximation and subsequent experience has shown that this method of 
generating overflow traffic is accurate and very useful in both simulations 
and analyses of traffic systems. 

I. INTRODUCTION 

A Poisson process which is alternately turned on for an exponentially 
distributed time and then turned off for another (independent) 
exponentially distributed time will be called an interrupted Poisson 
process-it can be viewed as a Poisson process modulated by a random 
switch. It was suggested by W. S. Hayward of Bell Laboratories that 
such a process be used to simulate overflow traffic. We will show that 
the interrupted Poisson process provides a simple and accurate method 
of simulating overflow traffic. 

The objective is to reduce the cost of computer simulations of 
traffic systems by using the interrupted Poisson process to model the 
overflow traffic. Generation of actual overflow traffic by simulating 
the behavior of the trunk group from which it overflows is time con- 
suming since a record must be kept of all calls which are offered to the 
subtending trunk group, whether they contribute to the overflow 
traffic or not. This is especially true when the traffic is overflowing a 
large trunk group. Moreover, the interrupted Poisson process provides 
a simple, approximate description of the overflow traffic and conse- 
quently facilitates analytical studies. 
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This method of generating overflow traffic has been used successfully 
in many studies of traffic systems. Examples of its application can be 
found in Refs. 1, 2, and 3 and also in numerous unpublished works. 
Since the method has been found to be very useful and requests for 
wider dissemination have been received by the author, this paper has 
been prepared. 

In Section II we derive the distribution of the number of busy 
trunks in an infinite trunk group when the offered traffic is generated 
by an interrupted Poisson process. The corresponding distribution for 
an overflow input has been computed by Kosten. 4 Now, matching the 
first three moments of the two distributions, we obtain equations for 
the parameters of the interrupted Poisson process. We also give these 
equations for a two-moment match. The errors committed in approxi- 
mating the distribution given by Kosten are also examined. In Section 
III we derive the interarrival time distribution for a traffic stream 
generated by an interrupted Poisson process. 

II. MOMENT-MATCH EQUATIONS 

Let the interrupted Poisson traffic be offered to an infinite trunk 
group. Let X be the intensity of the Poisson process, I/7 be the mean 
on-time of the random switch, 1/w be the mean off -time, and l/n be 
the mean service time. Let the state of the system be described by 
(m, n) with state probabilities p(m, n) where m is the number of 
servers busy, and n is the state of the switch taking on the value of 1 
or according to whether the process is on or off. 

The equilibrium equations for the stationary state probabilities are 

{my. + o>)p(m, 0) = yp(m, 1) + (jn + l)np(m + 1, 0), m ^ 0, 

(mn + 7 + X)p(w, 1) = cop(m, 0) + (m + l)np(m + 1, 1) 

+Xp(m- 1, 1), ro^l, (1) 

( 7 + X)p(0, 1) =o>p(0,0) +mp(1, 1). 

To solve this system of equations we introduce the probability generat- 
ing function 

G{z) = L p(m)z m = G^z) + G (z), 

where 

Gi(z) = L P(m, l)z m , G (z) = L p(m, 0)z m , 

and 

p(m) = p(m, 1) + p{m, 0) . 

Note that Gi{\) = Ep(m, 1) is simply the probability of the switch 
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being on and is given by a>/(y + a). Similarly, G (l) = E p(™, 0) 
= 7/(7 + <*>) is the probability of the switch being off. We will now 
derive the differential equations for Gi and G and obtain their 
solutions. 

From (1) and the definition of (7(2) we have 

M (z - l)G' (z) + «Go(2) - 7G1CO = 0, 
n(z - l)Gi(z) + (7 + X - Xz)Gi( Z ) - «G («) = 0. 

This system of equations is coupled. At the price of increasing the 
order, we can decouple the system by means of differentiation and 
simple substitution : 

M (2 - l)G'o' (2) + [ M + 7 + « - X(2 - 1)]G' («) 

- - uG (2) = 0, (2) 
M (2 - l)Gi' (2) + [ M 4- 7 + a; - X(2 - l)](?i(«) 

- - (a) + M )Gi(a) = 0. 
M 

Since we are interested in the moments of the distribution of the 
number of busy servers, it will be convenient to obtain solutions of (2) 
valid about 2=1. Subsequently, the series can be rearranged at the 
origin to yield the state probabilities. 

The change of variable 

X 
*--(«- 1) 

transforms (2) into 

Wo'(€) + (• - f)Q'o(f) - 0Qo(«) = 0. 



(3) 



where 



and 



Wi'tt) + (« - *)«!(«) - (i + /»)&(*) = 0. 



7 4- a; CO 

= 1 + , /? = -, 



Qi({) =G 1 (z), Q,(© =G 2 ( 2 ). 
Equations (3) have the solutions 

Qo(«) = CiFiOl; «; 0, 

O!(0 ^DAd + ^iie, 
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where 

- (a)„C 
ifi(a;6;c) = E — — : 

n =o(o) n n! 

with 

(a)„ = a(a + 1) (a 4- 2) • • • (a + n - 1), (a) = 1, 

is the confluent hypergeometric function and C and D are arbitrary 
constants. The conditions Gi (1) = 0/(7 + to) and Go(l) = 7/(7 + «) 
determine C and D and hence 



GOO - — — 1^1 

7 + O) 



> r \ 1 

(4) 



+ ifi 

7 + CO 



1 +/3;e;-(2 - 1) 



The power series in (z — 1) on the right-hand side of (4) is con- 
vergent since 1F1 is an entire function. The factorial moments of the 
distribution of the number of busy servers, being the coefficients of 
the expansion (4), are simply 

g(n) (1 ) =X " , \ , n = 0,1,2, ••-, (5) 

(7 + W)n 

where we have made the normalization n = I. For computational 
purposes, the following recurrence relation is useful : 

(to + n) 
G<»+i)(l) = X — '— (?(«>(1) . (6) 

(7 + « + n) 

To obtain the state probabilities p(m) = p(m,Q) + p(m, 1) we 
rearrange the series (4) at the origin and identify the coefficients in 
this new expansion as the state probabilities 

1 - (-1)* 

p(m) = — L <?<«+*>(l) 

wi!* =0 k 

.x; £(-»*- mi (7) 

w!>= m (j - to)! (7 + «)i 

We now show that the interrupted Poisson process provides an 
accurate method of generating overflow traffic. Let a erlangs of Poisson 
traffic be offered to an Erlang B system of c trunks. Let the overflow 
traffic be routed to an infinite trunk group and let Y be the number of 
busy servers in the infinite trunk group under statistical equilib- 
rium. The factorial moments 71/ (n) of Y and the state probabilities 
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f(m) = p[F = m"} have been computed by Kosten: 4 

tro(c) 



A/ (n) = a n 



a 



-(C)' 



/(m) = — £ 



(8) 



c\m\ k=ok\<j k+m (c)' 



where 

a c c /j + i—1 



<t (c)=-, c;(c) = Z( . )- -, i- 1,2,---. 

c ! , =0 \ i I (c — i) 

For a more accessible reference which gives the derivation of the 
factorial moments, see the appendix prepared by J. Riordan in Ref. 5. 

Now consider an interrupted Poisson traffic of original intensity X 
offered to an infinite trunk group, and let X be the number of busy 
servers under statistical equilibrium. The factorial moments of X and 
the state probabilities are given by (5) and (7) respectively. In this 
system we have three parameters, X, y, and a>, which are to be chosen 
so that the interrupted Poisson process gives the best approximation 
to the overflow process. In the present analysis we choose the moment 
approximation and, in particular, take the factorial moments. Thus we 
require that 

G<">(1) = Af (l0l n= 1,2,3, (9) 

and define the error 

E n {a, c) = f(?i) - p(n), n = 0, 1, 2, • • • . (10) 

With the aid of the recurrence relation (6), we can express (9) as 

to + n \ 

= a8 n , n = 0, 1, 2, (11) 

\u + 7 + n/ 
where 

M (n+1) a n (c) 

On = 



aM(n) <Tn+l(c) 

Equations (11) have the solution 

5 2 (5i — 8 ) — 5 (52 — 5i) 



X = a 



(Si - So) - (« 2 - 5i) 
So /X — a8i s 



u> = 



X \8i — So 



(12) 



o> /X — a8o\ 
" a \ So J 
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The parameters X, u, and 7 can be computed from (12) whenever 
a and c are specified. Often, however, one is concerned with final-route 
traffic in which only the mean, a, and variance, v, (or peakedness ratio 
z = v/a) are known. There are two ways to proceed in this case. 

First, using Wilkinson's equivalent random method, 5 determine S, 
the number of trunks, and A, the equivalent random load correspond- 
ing to the overflow traffic of mean a and variance v. Now set a = A 
and c = S and use eqs. (12) to compute X, a>, and 7 for a three-moment 
match. 

A second way to proceed is to determine the equivalent random 
load A as before, and set X = A in the last two equations of (12) to 
compute oj and y for a two-moment match. A satisfactory value of the 
equivalent random load A is given by Rapp's approximation : 6 

A = az + 3z(z - 1), (13) 

where z is the peakedness ratio v/a. In terms of a, z, and A, these 
equations can be written as 

a (A — a \ 

r .(i-i).. 

Note that for a two-moment match it is necessary to fix one of the 
parameters X, u, or 7. However, for a positive solution, they cannot 
be chosen arbitrarily-for positivity we must choose X > a, such as 
in (13). 

To illustrate the procedure, let us take an example with overflow 
traffic of mean 0.61 and variance 0.95. By Wilkinson's equivalent 
random method, we obtain S = 4 trunks and A = 3 erlangs. The 
first three moments from (8) with c = S = 4 and a = A = 3 are 

M a ) = 0.618 

M w = 0.708 (15) 

M w = 1.025. 

The three-moment match yields an interrupted Poisson process with 
parameters 

X = 2.553 

eo = 0.646 
7 = 2.022, 

and of course the same first three moments as in (15). 
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For the two-moment match, we set X = A = 3 and from (14) obtain 

co = 0.669 

7 = 2.621 . 

The first two moments will be equal to Af(u and Mm of (15) and the 
third moment is found to be 

G (3) (l) = 1.049 

which is not significantly different from M (3 ) in (15). Computing the 
state probabilities, we obtain the following result in which the two- 



moment match of the 


negative binomial fit 5 was 


included for 


comparison : 












Exact State 


Three-Moment Two-Moment 


Negative 


State 


Probabilities 


Match 


Match 


Binomial 





0.6164 


0.6161 


0.6149 


0.6086 


l 


0.2312 


0.2318 


0.2344 


0.2464 


2 


0.0965 


0.0962 


0.0953 


0.0924 


3 


0.0372 


0.0372 


0.0366 


0.0337 


4 


0.0130 


0.0130 


0.0129 


0.0122 


5 


0.0041 


0.0041 


0.0042 


0.0043 





0.0012 


0.0012 


0.0012 


0.0015 


7 


0.0003 


0.0003 


0.0003 


0.0005 


8 


0.0001 


0.0001 


0.0001 


0.0002 



We now examine the error E„(a,c). Since little additional com- 
putation is required to obtain the three-moment match, we feel it 
should be done to obtain a better fit. Consequently, we examined the 
error for a three-moment match only. Calculations of f(n), the state 
probabilities as they are given exactly, and p(n), the interrupted 
Poisson approximation, have been made. This was done for groups of 
1, 2, 4, 8, 16, 32, 40, and 48 trunks in the primary group with offered 
occupancies of 0.75 and 1, and for a group of 64 with offered occupancy 
of 0.75. For larger trunk groups, significant loss of accuracy prevented 
successful computation. We note that for the case of one trunk the 
approximation is exact. Where comparison could be made with pre- 
viously reported results for the negative binomial and the confluent 
hypergeometric approximations, 7 it was found that the three-moment 
match using the interrupted Poisson process gave uniformly better fit 
to the state probabilities. 

A typical result of the computations made is displayed in Table I. 
It was found that for a fixed occupancy the errors would increase, 
reach a maximum, and then decrease as the number of trunks was 
increased. This behavior can be seen in Fig. 1 where 

E = max \E„(a, c) \ = max | f(n) — p(n) \ 
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Table I — Three-Moment Match 
(c = 16, a/c = 0.75, X = 6.83, co = 0.366, 7 = 3.09) 









E n (12,16) 


E n (12,16) 


State n 


f(n) 


pin) 


= /(») - p(n) 


m 





0.64866359 


0.64765765 


0.00100593 


0.00155078 


1 


0.17173954 


0.17395867 


-0.00221912 


-0.01292144 


2 


0.08460030 


0.08381668 


0.00078362 


0.00926265 


3 


0.04523738 


0.04463748 


0.00059990 


0.01326127 


4 


0.02430285 


0.02417444 


0.00012841 


0.00528357 


5 


0.01280780 


0.01290307 


-0.00009527 


-0.00743819 


6 


0.00655960 


0.00668314 


-0.00012354 


-0.01883330 


7 


0.00325188 


0.00333032 


-0.00007845 


-0.02412386 


8 


0.00155792 


0.00158885 


-0.00003093 


-0.01985106 


9 


0.00072097 


0.00072378 


-0.00000280 


-0.00388859 


10 


0.00032234 


0.00031441 


0.00000793 


0.02460252 


11 


0.00013930 


0.00013020 


0.00000910 


0.06530564 


12 


0.00005822 


0.00005142 


0.00000681 


0.11694969 


13 


0.00002356 


0.00001937 


0.00000418 


0.17756000 


14 


0.00000923 


0.00000697 


0.00000226 


0.24472624 


15 


0.00000351 


0.00000240 


0.00000111 


0.31590026 


16 


0.00000129 


0.00000079 


0.00000050 


0.38876238 


17 


0.00000046 


0.00000025 


0.00000021 


0.46060149 


18 


0.00000016 


0.00000008 


0.00000009 


0.53013187 


19 


0.00000005 


0.00000002 


0.00000003 


0.59587179 



is plotted against c for different values of the offered occupancy a/c 
and again in Fig. 2 where the relative error corresponding to the state 
n at which \E n (a, c) | attained its maximum is plotted. 

III. INTERARRIVAL TIME DISTRIBUTION 

In analytical studies of systems, a description of overflow traffic is 
sometimes needed in terms of the interarrival time distribution. 2 This 
distribution is complicated and may be difficult to compute whenever 
the size of the trunk group which the Poisson traffic is overflowing is 
large. 8 If the interrupted Poisson traffic is used to generate the overflow 
traffic, the resulting interarrival time distribution, say A(t), is simple. 
Indeed, we show that this distribution is given simply by the mixture 
of two exponential distributions : 



A(t) = fci(l - «" ri + Ml - e- r, 0, 



(16) 



where 



n = §{X + « + 7 + V(X + co + T ) 2 -4M, 



r 2 = £{A + « + 7 - V(X + w + 7) 2 -4Xco}, 

X — r 2 

fci = , 

n - r 2 

fc 2 = l — M 
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10-5- 
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2 4 6 8 10 20 40 60 100 

SIZE OF PRIMARY TRUNK GROUP, C 

Fig. 1 — Maximum absolute error E. 



The proof goes as follows. Let W„ be the waiting time from t = 
until the time of the nth arrival and H n (t) be the distribution of W n . 
To obtain the interarrival time distribution, it is not necessary to 
find H„(t) for all re. The distribution Hi(t) and proper choice of initial 
conditions at t = is sufficient to find A (t) . We include the more 
general case here for completeness. 

If N(t) counts the number of arrivals in (0, t) and 



then 



Pk (t) = PIN® = fc], 



H n {t) = 1 - E P*(0 • 
*=0 



(17) 
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2 4 6 8 10 20 40 60 100 

SIZE OF TRUNK GROUP, C 



Fig. 2 — Relative absolute error corresponding to Fig. 1. 



This equation follows from the observation that W n > t if and only 
if N(t) < n. 

Taking the Laplace-Stieltjes transform of both sides of (17) we 
obtain 



where 



a n (s) = 1 - S £ TTk(s), 



a n (s) = / e-> l dH n {t), 



(18) 



Tr k (s) = / e- 8t p k (t)dt. 
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The initial conditions used to obtain (18) are 

1, k = 



P*(0+) = 



(19) 
10, fc > 0. 



We wall compute ir*(s) and hence determine a„(s). From the expression 
for a n (s) we then determine the interarrival distribution A (t). 

Let phm(t) be the probability that there were k arrivals in (0, t), 
given that an arrival occurred at t = and that at the instant t the 
switch is on if m = 1 and off if m = 0. These functions satisfy the 
system of differential equations 

poi(<) = «Poo(0 — (X + y)Poi(0j 

p'*i(0 = <*VM - (X + t)p*i(«) + Xp*_i.i(0, k = 1, 2, • • • , (20) 

p*o(0 = - wp* (<) + 7P*i(0» fc = 0, 1, 2, • • ■ , 

and the initial condition poi(0) = 1. 

Taking the Laplace transform of (20), we obtain 

S7Toi(s) = o>7r o(s) — (X + y)ir i(s) + 1, 

sxti(s) = wTfco(s) — (X + y)wki(s) + Xir*-i,i(s), fc = 1, 2, • • • , 

ST k0 (s) = — cowko(s) + Xtjh(s), fc = 0, 1, 2, • • • , 

where 

■"*»■(«) = / e-' l p k j(t)dt . 

This system of difference equations can be solved for 7r*o(s) and 
■n-ki(s) and hence 7r*(s), since 7rt(s) is the sum irko(s) + 7r*i(s). Omitting 
the details, we have 

s + co + 7 TX(s + u)l k 
»*(*) = /N /x , fc = 0, 1, 2, ■ • • , (21) 

t(«) L t(s) J 

where 

t(s) = s 2 + (X + 7 + cu)s + Xco. 

Substituting (21) into (18), we obtain 

TX(s + w)-|» , % 

«»(*) = /N , n = 1, 2, • • • . (22) 

L t(s) J 

The interarrival time distribution is given by the distribution of the 
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waiting time until the first arrival and hence its Laplace-Stieltjes 
transform is given by ai(s). Inverting «i(s), we obtain eq. (16). 
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